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The present status of the diquark model for exclusive reactions at moderately 
large momentum transfer is reviewed. That model is a variant of the Brodsky- 
Lepage approach in which diquarks are considered as quasi-elementary constituents 
of baryons. Recent applications of the diquark model, relevant to high energy 
physics with electromagnetic probes, are discussed: electromagnetic form factors 
of baryons in both the space-like and the time-like region, two-photon annihilations 
into proton-antiproton pairs as well as real and virtual Compton scattering. 

Exclusive processes at large momentum -transfer are described in terms of 
hard scatterings among quarks and gluons El. In this so-called hard scatter- 
ing approach (HSA) a hadronic amplitude is represented by a convolution of 
process independent distribution amplitudes (DA) with hard scattering ampli- 
tudes to be calculated within perturbative QCD. The DAs specify the distri- 
bution of the longitudinal momentum fractions the constituents carry. They 
represent Fock state wave functions integrated over transverse momenta. The 
convolution manifestly factorizes long (DAs) and short distance physics (hard 
scattering). The HSA has two characteristic properties, the power laws and 
the helicity sum rule. The first property says that, at large momentum transfer 
and large Mandelstam s, the fixed angle cross section of a reaction AB — > CD 
behaves as (apart from powers of log s) 



where n is the minimum number of external particles in the hard scattering 
amplitude. The power laws also apply to form factors: a baryon form factor 
behaves as 1/ Q"', a meson form factor as 1 /Q^. The counting rules are found to 
be in surprisingly good agreement with experimental data. Even at momentum 
transfers as low as 2 GeV the data seem to respect the counting rules. 
The second characteristic property of the HSA is the conservation of hadronic 
helicity. For a two-body process the helicity sum rule reads 



da/dt^f{e) s' 
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Figure 1: The Pauli form factor of the proton scaled by Q^. J^ata. are taken fromE. The 
sohd Hne represents the result obtained with the diquark modelLl. 

It appears as a consequence of utilizing the coUinear approximation and of 
dealing with (almost) massless quarks which conserve their helicities when in- 
teracting with gluons. The coUinear approximation implies that the relative 
orbital angular momentum between the constituents has a zero component in 
the direction of the parent hadron. Hence the helicities of the constituents sum 
up to the helicity of their parent hadron. The helicity sum rule is violated by 
20 — 30% in many cases. A particular striking example is the Pauli form factor 
of the proton which is determined by a helicity flip contribution. Its de- 
pendence (see Fig. |^) is compatible with a higher twist contribution (^ 
In explicit applications of the HSA (carried through only in valence quark ap- 
proximation and to lowest order QCD with very few exceptions) one encounters 
the difficulty that the data are available only at moderately large momentum 
transfer, a region in which non-perturbative dynamics may still play a crucial 
role. A general feature of these applications is the extreme sensitivity of the 
results to the DAs chosen for the involved hadrons. Only strongly end-point 
concentrated DAs provide results which are in a few cases in fair agreement 
with the dataB. These apparent successes of the HSA are only achieved at the 
expense of strong contributions from soft regions where one of the constituents 
carries only a tiny fraction of its parent hadron's momentum. This is a very 
problematical situation for a perturbative calculation. It should be stressed 
that none of the DAs used in actual applications leads to a successful descrip- 
tion of all large momentum transfer processes investigated so far. 
It seems clear from the above remarks that the HSA (in the valence quark 
approximation) although likely to be the correct asymptotic picture for exclu- 
sive reactions, needs modifications at moderately large momentum transfer. 
In a series of papers UlB such a modification has been proposed by us in 
which baryons are viewed as composite of quarks and diquarks. The latter are 
treated as quasi-elementary constituents which partly survive medium hard 
collisions. Diquarks are an effective description of correlations in the wave 
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functions and constitute a particular model for non-perturbative effects. The 
diquark model may be viewed as a variant of the HSA appropriate for moder- 
ately large momentum transfer and it is designed in such a way that it evolves 
into the standard pure quark HSA asymptotically. In so far the standard HSA 
and the diquark model do not oppose each other, they are not alternatives 
but rather complements. The existence of diquarks is a hypothesis. However, 
from experimental and theoretical approaches there have been many indica- 
tions suggesting the presence of diquarks. For instance, they were introduced 
in baryon spectroscopy, in nuclear physics, in astrophysics, in jet fragmentation 
and in weak interactions to explain the famous A/ =1/2 rule. Diquarks also 
provide a natural explanation of the equal slopes of jpeson and baryon Regge 
trajectories. For more details and for references, seeH. It is important to note 
that QCD provides some attraction between two quarks in a colour {3} state 
at short distances as is to be seen from the static reduction of the one-gluon 
exchange term. 

Even more important for our aim, diquarks have also been found to play a 
role in inclusive hard scattering reactions. The most obvious place to signal 
their presence is deep inelastic lepton-nucleon scattering. Indeed the higher 
twist contributions, convincingly observed by the NMClid, can be modelled as 
lepton-diquark elastic scattering. Baryon production in inclusive pp collisions 
o reveals the need for diquarks scattered clastically in the hard interaction 
For instance, kinematical dependences or the excess of the proton yield 
over the antiproton yield find simple explanations in the diquark model. No 
other explanation of these phenomena is known as yet. 

The diquark model: As in the standard HSA a helicity amplitude for the re- 
action AB CD is expressed as a convolution of DAs and hard scattering 
amplitudes (s, — t, —u ^ mf) 

M{s,t) ^ J dxcdxDdxAdxB^*ci^cWD{xD)TH{xi.,s,t)<^>AixA)'^B{xB) (3) 

where helicity labels are omitted for convenience. Implicitly it is assumed in 
that the valence Fock states consist of only two constituents, a quark and a 
diquark (antiquark) in the case of baryons (mesons). In so far the specification 
of the quark momentum fraction suffices; the diquark (antiquark) carries the 
momentum fraction I ~ Xi. If an external particle is point-like, e.g. a photon, 
the accompanying DA is to be replaced by (5(1 — Xi). As in the standard HSA 
contributions from higher Fock states are neglected. This is justified by the 
fact that that such contributions are suppressed by powers of ag/t as compared 
to that from the valence Fock state. 

In the diquark model spin (S) and spin 1 (V) colour antitriplet diquarks 
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are considered. Within flavour SU(3) the S diquark forms an antitriplet, the 
V diquark an sixtet. Assuming zero relative orbital angular momentum be- 
tween quark and diquark and taking advantage of the coUinear approximation, 
the valence Fock state of a ground state octet baryon B with helicity A and 
momentum p can be written in a covariant fashion (omitting colour indices) 

\B;p,X) = fs'P^{x)Bsuip,X) + fv^^{x)Bvh''+p"/mBhrouip,X)/V3 (4) 

where u is the baryon's spinor. The two terms in (Q) represent configurations 
consisting of a quark and either a scalar or a vector diquark, respectively. The 
couplings of the diquarks with the quarks in a baryon lead to flavour functions 
which e.g. for the proton read 

Ss = w5'[„,d] By ^[uV^uM}-V2dViu,u}]/V3. (5) 

The DAs ^s{v) conventionally normalized as / da;$ = 1. The constants fs 
and fv play the role of the configuration space wave functions at the origin. 
The DAs containing the complicated non-perturbative bound state physics, 
cannot be calculated from QCD at present. It is still necessary to parameterize 
the DAs and to fit the eventual free parameters to experimental data. Hence, 
both the models, the standard HSA as well as the diquark model, only get 
a predictive power when a number of reactions involving the same hadrons 
is investigated. In the diquark model the fol|k)"[5ii|ng DAs have been proven to 
work satisfactorily well in many applications QtJ: 

$f (x) =iV|x(l - exp [-b'{ml/x + m|/(l - x))] (6) 
^y{x)=N^x{l - x)\l + b.8x- I2.bx^) exp [-^^(m^/a; + m^/(l - x))] . 

The constants iVj' and Ny are fixed through the normalization (e.g. for the 
proton = 25.97, = 22.92). The DAs exhibit a mild flavour dependence 
via the exponential whose other purpose is to guarantee a strong suppression 
of the end-point regions. The masses in (^) are constituent masses; for u and 
d quarks we take 350 MeV and for the diquarks 580 MeV. Strange quarks and 
diquarks are assumed to be 150 MeV heavier that the non-strange ones. It is 
to be stressed that the quark and diquark masses only appear in the DAs (^; 
in the hard scattering kinematics they are neglected. The final results (form 
factors, amplitudes) depend on the actual mass values mildly. The transverse 
size parameter b is fixed from the assumption of a Gaussian transverse mo- 
mentum dependence of the full wave function and the requirement of a value 
of 600 MeV for the mean transverse momentum (actually b — 0.498 GeV~^). 
As the constituent masses the transverse size parameter is not considered as a 



4 



free parameter since the final results only depend on it weakly. 
The hard scattering amplitudes Tu, determined by short-distance physics, are 
calculated from a set of Feyman graphs relevant to a given process. Diquark- 
gluon and diquark-photon vertices appear in these graphs which, following 
standard prescriptions, are defined as 

SgS: igst" {,Pi+P2)^. 

VgV : -i gst" {jjapipi + P2)t, - ff/?^ [(1 + «) P2 - 

[(1 + - KP2]^ } (7) 

where Qs = \f^Kcxl is the QCD coupling constant, k is the anomalous magnetic 
moment of the vector diquark and t'^ = A°/2 the Gell-Mann colour matrix. For 
the coupling of photons to diquarks one has to replace gst"^ by —\J AnaeB where 
a is the fine structure constant and eu is the electrical charge of the diquark 
in units of the elementary charge. The couplings DgD are supplemented by 
appropriate contact terms required by gauge invariance. 

The composite nature of the diquarks is taken into account by phenomenolog- 
ical vertex functions. Advice for the parameterization of the 3-point functions 
(diquark form factors) is obtained from the requirement that asymptotically 
the diquark model evolves into the standard HSA. Interpolating smoothly be- 
tween the required asymptotic behaviour and the conventional value of 1 at 
= 0, the diquark form factors are actually parametrized as 



Pf^^') = 7^.^ ^rW^)- 7,2^) • (8) 



>s^^- \Qv + Q" 

The asymptotic behaviour of the diquark form factors and ±h.e connection to 
the hard scattering model is discussed in more detail in Ref. ErQ. In accordance 
with the required asymptotic behaviour the n-point functions for n > 4 are 
parametrized as 

Fi")(g2) = a5Ff (g^), ^v^W') = («^Qr^) 4'\Q')- (9) 

The constants as,v are strength parameters. Indeed, since the diquarks in 
intermediate states are rather far off-shell one has to consider the possibility 
of diquark excitation and break-up. Both these possibilities would likely lead 
to inelastic reactions. Therefore, we have not to consider these possibilities ex- 
plicitly in our approach but excitation and break-up lead to a certain amount 
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of absorption which is taken into account by the strength parameters. Admit- 
tedly, that recipe is a rather crude approximation for n > 4. Since in most 
cases the contributions from the n-point functions for n > 4 only provide small 
corrections to the final results that recipe is sufficiently accurate. 
Special features of the diquark model: The diquark hypothesis has striking 
consequences. It reduces the effective number of constituents inside baryons 
and, hence, alters the power laws (|l|) . In elastic baryon-baryon scattering, for 
instance, the usual power s^^*^ becomes s~^F(s) where F represents the net 
effect of diquark form factors. Asymptotically F provides the missing four 
powers of s. In the kinematical region in which the diquark model can be 
applied {—t, —u > 4GeV^), the diquark form factors are already active, i.e. 
they supply a substantial s dependence and, hence, the effective power of s 
lies somewhere between 6 and 10. The hadronic helicity is not conserved in 
the diquark model at finite momentum transfer since vector diquarks can flip 
their helicities when interacting with giuons. Thus, in contrast to the standard 
HSA, spin-flip dependent quantities like the Pauli form factor of the nucleon 
can be calculated. 

Electromagnetic nucleon form factors: This is the simplest application of the 
diquark model and the most obvious place to fix the various parameters of the 
model. The Dirac and Pauli form factors of the nucleon are evaluated from 
the convolution formula (||) with the DAs (^) and the parameters are deter- 
mined from a best fit to the data in the space-like region. The following set of 
parameters 

fs = 73.85 MeV, Q| = 3.22 GeV^ as = 0.15, , . 

/y = 127.7 MeV, Q^ = 1.50GeV^ ay = 0.05, k = 1.39 ; ^ ' 

provides a good fit of the dataB. as is evaluated with Aqcd = 200 MeV and 
restricted to be smaller than 0.5. The parameters Qs and Qy, controlling the 
size of the diquarks, are in agreement with the higher-twist effects-observed 
in the structure functions of deep inelastic lepton-hadron scattering EJ if these 
effects are modelled as lepton-diquark elastic scattering. The Dirac form factor 
of the proton is perfectly reproduced. The results for the Pauli form factor are 
shown in Fig. 0. The predictions for the two neutron form factors are also in 
agreement with the data. However, more accurate neutron data are needed in 
the region of interest in order to determine the model parameters better. 
The nucleon's axial form factor lI and its electromagnetic form factors in the 
time- like regions □ have also been evaluated. Both the results compare well 
data. Even electroexcitation of nucleon resonances has been investigated 
. In the-case of the A^A form factor the model results agree very well with 
recent datata. 
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Figure 2: Typical Feynman graphs contributing to 7(*'p 7p. 
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Figure 3: The scaled cross section for RCS off p«otons vs. cos 8 for three different photon 
energies. The experimental data arc taken fromF^I. 



Real Compton scattering (RCS): 7p — s- 7p is the next reaction to which the 
diquark model is apphed to. Since the only hadrons involved are again protons 
RCS can be predicted in the diquark model without any adjustable parameter. 
Typical Feynman graphs contributing to that process are shown in Fig. |^. The 
results of the diaiiark model for RCS are shown in Fig. || for three different 
photon energies EIEI. Note that in the very forward and backward regions the 
transverse momentum of the outgoing photon is small and, hence, the diquark 
model which is based on perturbative QCD, is not applicable. Despite the 
rather small energies at which dataEj are available, the diquark model is seen 
to work rather well. The predicted cross section does not strictly scale with 
s~^. The results obtained within the standard HSA are of similar qualitylHl. 
The diquark model also predicts interestiiig photon asymmetries and spin cor- 
relation parameters (see the discussion inQ). Even a polarization of the proton, 
of the order of 10%, is obtained H. This comes about as a consequence of hc- 
licity flips generated by vector diquarks and of perturbative phases produced 
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Figure 4: The integrated my pp cross section (| cos 6 1^-0.6). The sohd hne represents the 
diquark model prediction cl Data are taken from CLEOlij. 

by propagator poles appearing within the domains of the momentum fraction 
integrations. The poles are handled in the usual way by the le presription. 
The appearance of imaginary Jiarts to leading order of as is a non-trivial pre- 
diction of perturbative QCDE3; it is characteristic of the HSA and is not a 
consequence of the diquark hypothesis. 

Two-photon annihilation into pp pairs: This process is related to RCS by cross- 
ing, i.e. the same set of Feynman graphs contributes (see Fig. ||). The only 
difference is that now the diquark form factors are needed in the time-like 
region. The expressions (^,^ represent an effective parameterization of them 
valid at space- like Q^- Since the exact dynamics of the diquark system is not 
known it is not possible to continue these parameterizations to tke time-like 
region in a unique way. A continuation can be defined as follows is re- 

placed by —s in (H,^ guaranteeing the correct asymptotic behaviour and, in 
order to avoid the appearance of unphysical poles at low Q^, the diquark form 
factors are kept constant once their absolute values have reached cq ~ 1.3 H. 
The same definition of the time- like diquark form factors is used in the analysis 
the proton form factor in the time-like region. The diquark model predictions 
for the integrated 77 pp cross section is compared to the CLEO datall^ in 
Fig. ^. At large energies the agreement between predictions and experiment is 
good. The predictions for the angular distributions are in agreement with the 
CLEO data too. The standard HSA on the other handairedicts a cross section 
which lies about an order of magnitude below the dataBy. Recently CLEO has 
also measured two-photon annihilations into AA pairs til. At large energies the 
integrated AA cross section is about a factor of 2 smaller than that one for 
annihilations into pp pairs. This may be taken as a hint at a more complicated 
A wave function than (||j6|). 
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Virtual Compton scattering (VCS): This process is accessible through ep ~* 
epj. An interesting element in that reaction is that, besides VCS, there is 
also a contribution from the Bethe-Heitler (BH) process where the final state 
photon is emitted from the electron. Electroproduction of photons offers many 
possibilities to test details of the dynamics: One may measure the s, t and 
dependence as well as that on the angle between the hadronic and leptonic 
scattering planes. This allows to isolate cross sections for longitudinal and 
transverse virtual photons. One may also use polarized beams and targets and 
last but not least one may measure the interference between the BH and the 
VC contributions. The interference is sensitive to phase differences. 
At s, —t and —u ^ (or small | cos 6*1 where 6 is the scattering angle of 
the outgoing photon in the photon-proton center of mass frame) the diquark 
model can also be applied to VCS u. Again there is no free parameter in 
that calculation. The relevant Feynman graphs are the same as for RCS (see 
Fig. ^). The model can safely be applied for s > 10 GeV^ and | cos6'| < 0.6. 
For the future CEBAF beam energy of 6 GeV the model is at its limits of 
applicability. However, since the diquark model predictions for real Compton 
scattering do rather well agree with the data even at s > 5 GeV^ (see Fig. ||) 
one may expect similarlv^good agreement for VCS. Predictions for the VCS 
cross section are given ina. The transverse cross section (which, at — 0, is 
the cross section for RCS) is the dominant piece. The other cross sections only 
become sizeable for larger values of | cos 6*1. Examination of the Bethe-Heitler 
contribution to the process ep — *■ epj reveals that it is small as compared to 
the VCS contribution at high energies, small values of | cos 6*1 and for an out- 
of-plane experiment, i.e. cj) > 50°. 

Of interest is also the electron asymmetry in ep —^ epj: 

A, = 'Jil-fid (11) 

cr(+) + (t(+) 

where ± indicates the helicity of the incoming electron. Al measures the 
imaginary part of the longitudinal - transverse interference. The longitudinal 
amplitudes for VCS turn out to be small in the diquark model (hence AY^'-' is 
small). However, according to the model, Al is large in the region of strong 
BH contamination (see Fig. |^). In that region, A^ measures the relative phase 
(being ofjjerturbative origin from on-shell going internal gluons, quarks and 
diquarksEa) between the BH amplitudes and the VCS ones. The magnitude 
of the effect shown in Fig. |^ is sensitive to details of the model and, therefore, 
should not be taken literally. Despite of this our results may be taken as an 
example of what may happen. The measurement of Al, e.g. at CEBAF, will 
elucidate the underlying dynamics of VCS strikingly. 



9 



s=5GeV^ Q'=1GeV^ k,L=6GeV 

0.50 
0.30 
0-10 



-0.10 

20 40 60 80 100 120 140 160 180 
(degres) 

Figure 5: Diquark model predictions for the electron asymmetry in ep — > 6^7!. 

Summary and outlook: The diquark model which represents a variant of the 
HSA, combines perturbative QCD with non-perturbative elements. The di- 
quarks represent quark-quark correlations in baryon wave functions which 
are modelled as quasi-elementary constituents. This model has been applied 
to many photon induced exclusive processes at moderarely large momentum 
transfer (typically > 4GeV^). From the analysis of the nucleon form fac- 
tors the parameters specifying the diquark and the DAs, are fixed. Compton 
scattering and two-photon annihilations of can then be predicted. The com- 
parison with existing data reveals that the diquark model works quite well and 
in fact much better then the pure quark HSA. 

Predictions for the VCS cross section and for the ep ejry cross section have 
also been made for kinematical situations accessible at the upgraded CEBAF 
and perhaps at future high energy accelerators like ELFE@HERA. According 
to the diquark model the BH contamination of the photon electroproduction 
becomes sizeable for small azimuthal angles. The BH contribution also offers 
the interesting possibility of measuring the relative phases between the VC 
and the BH amplitudes. The phases of the VC amplitudes are a non-trivial 
phenomenon generated by the fact that some of the internal quarks, diquarks 
and gluons may go on mass shell. The electron asymmetry is particularly 
sensitive to relative phases. In contrast to the standard HSA the diquark model 
allows to calculate helicity flip amplitudes, the helicity sum rule (||) does not 
hold at finite . One example of an observable controlled by helicity flip con- 
tributions is the Pauli form factor of the proton. Also in this case the diquark 
model accounts for the data. 

Photo- and electroproduction of mesons can also be calculated within the di- 
quark model. However, these reactions are already quite complicated since all 
together 158 Feynman graphs contribute. The analysis of this class of pro- 
cesses is not yet finished. First results exist only for 7p — > -fr'^*'At2l. 
The only discrepancy between data and diquark model predictions known to 
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us, is to be seen in the process 77 AA. This failure seems to indicate that the 
A wave function, and perhaps those of other hyperons are more comphcated 
than (0,^). A solution of this problem may necessitate a combined analyis of 
77 AA and K'-*'>A. 
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